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The c-axis reflectivity spectrum of underdoped YBa2Cu306.6 (YBCO) is measured below T c = 
59K in parallel magnetic fields _ff//Cu02 up to 7 T. Upon application of a parallel field, a new 
peak appears at finite frequency in the optical conductivity at the expense of suppression of c- 
axis condensate weight. We conclude that the dramatic change originates from different Josephson 
coupling strengths between bilayers with and without Josephson vortices. We find that the 400cm -1 
broad conductivity peak in YBCO gains the spectral weight under parallel magnetic field; this 
indicates that the condensate weight at ui — is distributed to the intra-bilayer mode as well as to 
the new optical Josephson mode. 

PACS numbers: 74.25.Gz, 74.72. Bk, 74.50.+r 



Josephson plasma oscillation, which generally appears 
in the E//c optical response of high T c cuprates is a 
direct consequence of the two-dimensionality of their su- 
perconductivity |il 0, p|, El . Depending on the anisotropy 



parameter 7 



Xr.rx. 



the Josephson plasma fre- 



quency ranges from the microwave (7 ~ 300 — 1000 
for Bi 2 Sr 2 CaCu20 8+( 5 (Bi-2212) §, |, §) to the far 
infrared frequency (7 ~ 50 for underdoped YBCO 
§1,0,0,12) and La 2 _ x Sr x Cu0 4 (LSCO) The 
optical/microwave response originates from the restora- 
tion of coherent transport due to the tunneling of su- 
perconducting carriers between the Cu02 layers. The 
tunneling current, which is driven by the phase-shift be- 
tween the CuC>2 layers, is most directly perturbed by the 
magnetic field applied parallel to the layers; the inter- 
play between parallel field and Josephson tunneling has 
recently become an interesting subject. Since the tun- 
neling current is involved in the screening of the parallel 
field, vortices have a large size of ~ 7s within the lay- 
ers, whereas in the c-axis direction, they are squeezed 
to the separation s of the superconducting layers. This 
very anisotropic vortex in the parallel field is called a 
Josephson vortex. 

In Bi-2212, there have been microwave resonance mea- 
surements under parallel applied magnetic fields (l4|, [II], 
|l§, and resonance peaks have been observed in mi- 
crowave absorption spectra. However, in the microwave 
study, it is not possible to discuss the origin of the ob- 
served modes quantitatively. In this Letter, we report a 
dramatic change of the Josephson plasma mode and an 
emergence of new modes under parallel fields in under- 
doped YBCO, where relatively small anisotropy 7 ~ 50 
brings Josephson resonance frequencies up to the far in- 
frared regime. The technique of optical reflectivity was 
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employed which enabled us to investigate both the reso- 
nance energies of longitudinal and transverse modes and 
their oscillator strengths, necessary for a quantitative dis- 
cussion of the sum rule. 

Another key issue of the c-axis response of YBCO is the 
effect of the bi-layer structure. Observations of a broad 
peak at ~ 400cm -1 in the optical conductivity which ap- 
pears below the pseudo-gap temperature have been made 
by several groups §, §, 0, ||, ||. This "400cm -1 " 
peak has been interpreted as the out-of-phase coupling 
of the two types of Josephson junctions in YBCO, which 
correspond to the intra- and the inter-bilayer coupling 
Jl7| , Eq . Similar observations of this "2nd Josephson 
plasma mode" have been made in Bi-2212 with the 
same double pyramid structure as YBCO and T*-type 
compounds, such as SmLai_ x Sr x Cu04 [|l9| ^0], ^J, which 
also has two types of insulating layers. Since Josephson 
vortices penetrate more easily into junctions with weaker 
Josephson coupling, the applied field is expected to pri- 
marily modulate the inter-bilayer coupling in YBCO. It 
is also interesting to explore the effect of parallel field on 
the 400cm -1 mode. If there were any effect, it would be 
direct evidence that superconductivity is involved in the 
formation of this mode. 

The specimen we measured is a YBa2Cii306.6 single 
crystal grown by the pulling method Q . A large crystal 
(10 x 10 x 10 mm 3 ) was cut along the a x c-plane, post 
annealed in latm O2 gas at 750°C in order to ensure 
the oxygen content, and then quenched to the room tem- 
perature. Magnetic susceptibility (not shown) exhibits a 
sharp superconducting transition at T c = 59K, with the 
transition width of less than 2K. The crystal was pol- 
ished using AI2O3 powders, and varnished on a sample 
holder with two 06 mm holes, together with a gold coated 
glass mirror. We employed a Fourier-transform infrared 
spectrometer IFS-113v (Bruker) which is equipped with a 
superconducting magnet SpectroMag (Oxford). We per- 
formed calibration measurements in a regular cryostat in 
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FIG. 1: c-axis reflectivity of YBCO, under an applied mag- 
netic field H I /afe-plane. All the data shown are taken after 
zero-field cooling, except the data labeled 7T FC (field cool- 
ing)- 



zero-field. The result agreed with that measured in the 
magnet. Small misalignments of the magnetic field with 
respect to the CuC>2 plane direction do not have serious 
effects on the result, as the anisotropy is much smaller 
than that of Bi-2212. 

In Figsj^(a) and [l](b), we show the c-axis optical re- 
flectivity spectra of YBCO at 5.5K and various fields. 
At zero magnetic field, the Josephson plasma edge is lo- 
cated at ~ 45cm" 1 and a broad bump at ~ 400cm is 



observed as has been previously reported || ||. [n], 11,12). 
When a parallel magnetic field {H a b > 1 Tesla) is applied 
after the sample has been cooled to 5.5K in zero-field 
(zero-field cooling: ZFC), the reflectivity edge shifts to 
higher energy, but simultaneously the reflectivity in the 
low energy side rapidly decreases, forming a peak. This 
result demonstrates a surprising sensitivity to the paral- 
lel magnetic field. By contrast, a slight and monotonic 
low energy shift is observed for fields applied perpendicu- 
lar to the planes (not shown). Fig.jl] also shows the data 
taken after field-cooling at 7T (labeled as "7.0T FC"). 
In this field-cooling situation, the reflectivity increases 
below ~ 25cm" 1 forming an additional edge, and the lo- 
cation of the higher energy peak and edge structure is 
different from that at 7T(ZFC). The spectral difference 
between FC and ZFC indicates that the observed spectral 
change is related to the pinning of Josephson vortices. 

To make a quantitative analysis of the data, we em- 
ployed the Kramers-Kronig analysis of the reflectiv- 
ity, and calculated the optical conductivity a(uj) — 
u>Im(e) /Air and the energy loss function Im(—l/e). The 
results are shown in Figs.Q(a)-(c). A remarkable change 
in the spectrum is the appearance of a peak in the optical 
conductivity, labeled as "cr" in Fig.^(a), which shifts to 
higher energies with increasing field. The peak intensity, 
on the other hand, does not depend much on the field 
strength. The cr-mode obtains appreciable intensity al- 
ready at H ao = 2T (probably at an even lower field in 



30 L 



50 

~ 40 
> E 
o 

§ 35 

~° X 20 

° w -m 
O b 10 





c 
o 

3jT 0.1 

1 M0.05 
o 



YBaCuO 5.5KH//ab 

2 3 6.6 







c 
o 

o 5? 0.1 

% - 0.05 
o 



70 

" 60 
E 

3 50 
I 40 
| 30 
% 20 
I 10 




L (a) ZFC 5.5K 


- 


7T J 

: 

0T : 

i , , , : 


! (b) ZFC 5.5K p 

; I// 

: , , ,dE=r*^ 


- i i i | i i i | i i i | i i i 

(c)FC5.5K p 

: i i i 1— ri 


7T : 

5T \ 

\ 

0T j 

i i i i : 



20 



40 



60 



Wave number (cm 1 ) 



80 



100 



:(d)ZFC5.5K ] 

:|m(-1/e) — . I 


: (e) FC 5.5K * : 

Mm(^1/£)__^ " i 


^ ' " (B=(a nVJH - 

: ^ JO: 
^ H o =<J> o /7S 2 (-^50); 


I ° ^^"^^ 



01234567 01234567 

Parallel magnetic field H ab (T) 



FIG. 2: Optical conductivity (a) and energy loss function (b) 
at 5.5K after zero-field cooling, (c) Energy loss function at 
5.5K after field cooling. Measured peak position of (d) ZFC 
and (e) FC. The dashed line in (d) and (e) is a theoretical 
prediction for the high-field limit based on Ref. |27]]. 



view of the data at IT in Fig.|](a)), which suggests an 
abrupt growth of intensity upon application of a parallel 
field. In FigSJ^(b) and |^(c), the energy-loss functions 
which generally measure longitudinal modes are com- 
pared for zero- field cooling (ZFC) and field-cooling (FC). 
In ZFC [Fig|](b)], there is only one peak observed in our 
experimental window (to > 15cm -1 ). In FC [FigJ|(c)], 
there are two peaks emerge within our window. The new 
longitudinal mode (labeled "a") increases its intensity 
with H a i,. The original longitudinal mode (labeled "/?", 
corresponding to the Josephson plasma mode at H ao = 0) 
decreases in intensity, showing that the spectral weight 
transferred from j3- to a-mode. This happens also for 
ZFC, though the a-peak is not seen above 15cm _1 . The 
energy positions of the peaks are plotted in Figs]|(d) and 
0(e) as a function of H ao . The peak energy of the /3-mode 
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FIG. 3: (a) c-axis optical conductivity of YBCO, under 
a parallel field H a b. The rectangles indicate the weight of 
low-energy condensate as derived from fitting to the function 
Re(e) « Eoo - p s /u 2 shown in the inset for 0T and 7T(FC) 
at 5.5K. The filled circles in the inset indicate the data points 
used for the fitting, (b) Optical conductivity on a wider scale. 



shifts upward, whereas the peak energy of the a-mode 
decreases slightly with H a t- 

Existence of longitudinal a-mode, or equivalently, the 
reflectivity edge below the peak energy, indicates that 
there is another transverse mode at ui a > < 15cm -1 . This 
cr'-mode may be either an u> = superconducting conden- 
sate or a dissipative mode at finite energy. We estimated 
the spectral weight of this cr'-mode (p s ) by analyzing the 
low-energy asymptotic behavior of the dielectric function 
Re(e) w Eoo — Ps/lo 2 (see inset of Fig.||(a)). This analysis 
is exact for superconducting condensate, but may have 
~ ±30% error in the spectral weight if uv is finite. For 
H a b = p s is a genuine c-axis condensate weight, and es- 
timated to be p s = 4.0(1) xlO 4 cm" 2 . For H ab = 7T(FC), 
there remains a weight of 1-5(4) x 10 4 cm -2 at low energy 
which locates either at w = or at a finite energy uj a i. 
The low energy condensate weights are shown in Figj^(a) 
as the height of the rectangles with a width of 10cm -1 . 
For measurements at 7T(ZFC), 1/oj 2 behavior in Re(e) 
is no longer seen, which suggests that the low-w weight 
is very small. 

Optical conductivity over wider energy range is also 
shown in Fig.||(b). It is clear that the conductivity peak 
at 400cm _1 gains intensity under magnetic fields, espe- 
cially at 7T(ZFC). The integrated intensity of the cr-peak 
and the increase in the 400cm~ 1 peak intensity yield 
S(u} a ) = 1.8(1) x 10 4 cm" 2 and AS(u T ) = 2.2(2) x 10 4 
cm -2 respectively at 7T(ZFC), which completely ex- 
hausts the weight of superconducting condensate at zero- 
field. On the other hand, at 7T(FC), S(u> a ) = 1.1(1) x 
10 4 ciTr 2 and AS(lj t ) = 0.5(2) x 10 4 cm~ 2 . The sum ex- 
plains 60-70% of the missing p s (~ 2.5(4) x lO^mT 1 ) 
transferred to the cr-mode and the 400cm- 1 peak. Sub- 
stantial weight remains in the cr'-mode, which is also evi- 
denced by the evolution of the ct-peak in Im{—l/e) with 



increasing H a t seen in Fig. ||(c). The sum rule analy- 
sis gives evidence for spectral-weight transfer from the 
superconducting condensate to the other two or three 
finite-u; transverse modes, enforced by parallel fields or 
by modulation of the inter-/intra-bilayer coupling. 

Now we discuss the origin of the cr-mode. It is proposed 
that Josephson vortex lattice has a quantized c-axis lat- 
tice parameter (= ./V x c), whereas its in-plane lattice 
parameter is tunable to match the external field p5[ . 
In this situation, the inter-bilayer coupling is expected 
to be modulated, differentiating the Josephson coupling 
strength between bilayers with and without Josephson 
vortices, giving rise to a transverse (optical) Josephson 
plasma mode at finite u>. For underdoped YBCO, the 
anisotropy parameter 7 ~ 50 scales the characteristic 
field H = $o/7s 2 to « 30T ^ §§. Calculations of free- 
energies show that Josephson vortices exist in one out of 
every two or more layers when H ab < H* = Ho/3 « 10T 
p5| , [26| , which very nearly matches our experimental con- 
ditions. A phenomenological model which describes the 
optical response of a Josephson junction array with more 
than two kinds of insulating layers has been proposed by 
van der Marel et al. which is called as "multi-layer 
model" . This model has been employed for the anal- 
ysis of c-axis optical resp onse of bi-layer compounds in 
zero-field (YBCO @ @ and Bi-2212 j§3) and T* mate- 
rials j^, |2(], [2l[] where two kinds of Josephson couplings 
are expected due to the crystal structure. We model 
YBCO in parallel magnetic fields by a simplified multi- 
layer model which takes into account only the low energy 
modes in the energy range (u> < 80cm -1 ) where optical 
phonons are absent. We assume that the possible cr'- 
mode is located at uj = 0: 



e(w) (^ 2 - ^)(^ 2 



UJ 2 (liJ 2 — LO 2 ) 



(1) 



After a simple calculation, the spectral weight of the 
transverse modes S(ui a ) and p s are given by a function 
of the energies of these resonance modes. At H a {, = 
7T(FC), the estimate is straightforward, since all the 
relevant resonance frequencies are observed. They are 



23.6cm 



44.4cm 1 and ujp = 54.0cm 



Using these frequencies and the dielectric constant = 
21, the calculation gives S(w ff ) = 1.4 x 10 4 cm~ 2 and 
p s = 1.7 x 10 4 cm- 2 , in reasonable agreement with the 
observed cr-peak intensity S(oj a ) — 1.1(1) x 10 4 cm -2 and 
the cr'-peak intensity (or superconducting condensate) 
p s = 1.5(4) x 10 4 cm" 2 . At 7T(ZFC), the location of 
the a-peak is uncertain which makes the estimate some- 
what ambiguous. Also in this case, the observed cr-peak 
intensity 1.8(1) x 10 4 cm -2 is consistent with a rough es- 
timate of w 2 x 10 4 cm- 2 . These agreements supports the 
hypothesis that the cr-mode originates from the periodic 
modulation of inter-bilayer Josephson coupling strength 
due to the Josephson vortices. 

A theoretical investigation of Josephson plasma in par- 
allel field was performed by Bulaevskii et al. under 
the ideal condition of absence of pinning to the Joseph- 
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son vortices |27|] , They found two strong resonance 
modes: one corresponding to "vortex sliding mode" , and 
the other arising from the Josephson plasma dispersion 
folded to the zone center by the formation of a vor- 
tex superlattice ^7|. The spectral weight is transferred 
from the latter to the former with H a b- In real materi- 
als, the vortex sliding mode is expected to be influenced 
by pinning of the Josephson vortices; the tr'-mode (or 
the associated longitudinal a-mode) in YBCO depends 
on the field-cycle and may be assigned to the pinned 
sliding mode of Josephson vortices. The energy of the 
"folded Josephson plasma mode" calculated in high field 
limit Jp7|, presents a linear field dependence as shown in 
Figs.g(d) and ||(e) by the dashed line. It appears that 
the energy of the cr-mode (or the associated /3-mode) in 
YBCO tends to approach this line at higher fields. In 
fact, the spectral weight transfer between a- and /3-mode 
shown in Figs. 0(b) and ^(c) is consistent with the the- 
oretical result p7| . Recent microwave measurement of 
Bi-2212 under parallel field has revealed two resonance 
modes Jil|. In the light of the present results for YBCO, 
it is likely that the higher energy mode corresponds to 
the /?- (or a-) mode, and the low-energy mode to the 
a-mode (or cr'-mode). 

The reason why the transverse a- and 400cm~ 1 modes 
are considerably more enhanced in ZFC than in FC might 
be related to the pinning of Josephson vortices. As the 
alternating stack of insulating layers with and without 
Josephson vortices is most likely the source of the cr- 
mode, disorder in stacking would reduce the cr-mode in- 
tensity. This seems to be the case with the cr-mode for 
FC, since random distribution of pinning centers tends to 



make the stacking of Josephson vortices more disordered. 
In a T* cuprate Nd2-x-j/Sr x Cej / Cu04, reduced intensity 
of the transverse Josephson plasma mode has been inter- 
preted as the result of disorder in the crystalline structure 

In LSCO, no new transverse cr-mode under parallel 
field is seen fl2S|, g |o) 



although the anisotropy parame- 
ter 7 is in the same order of magnitude with underdoped 
YBCO. Only a low-w shift of the Josephson plasma edge 
is observed with increasing parallel field, suggesting that 
the Josephson vortices would be highly disordered seem- 
ingly due to fluctuations of stripe order in LSCO. 

In conclusion, we have demonstrated a remarkable 
effect of parallel magnetic field on the c-axis Joseph- 
son plasma mode in underdoped YBCO. Parallel mag- 
netic field modulates the inter-bilayer Josephson coupling 
strength along the c-axis due to the presence of inequiv- 
alent insulating layers with and without Josephson vor- 
tices. As a result, one optical (transverse) mode appears, 
corresponding to the anti-phase Josephson current oscil- 
lations between two inequivalent junctions. Another con- 
sequence of the modulation of the inter-bilayer coupling 
and/or suppression of inter-bilayer coherent pair hopping 
is the enhanced intensity of the 400cm _1 mode. This 
gives evidence that this mode is a marker of supercon- 
ducting Josephson coupling within a bilayer. 
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